NASA TECHNICAL MEMORANDUM 


/N 3 ? 

NASA TM 88491 


ON CONTACT PROBLEMS OF ELASTICITY THEORY 
A. I. Kalandiya 



Translation of "K kontaktnym zadachim teorii uprugosti," IN 
Prikladnaya mateniatika i mekhanika, Vol. 21, No. 3, 1957, 
Moscow, pp. 389-398. 


(NASA -Til-8849 1) CN CCNTAC1 EICBIEMS OF N37-14729 

ELASTICITY THEOFY (National Aeronautics and 
Space Administration) 20 p CSCL 46E 

Unclas 

G3/39 43669 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 OCTOBER, 1986 



SI »HDtRD llH[ Pag [ 


^ po« i Nt. 

NASA TM— 88491 


^ • Till* «rt<2 



ON CONTACT PROBLEMS OF ELASTICITY 
THEORY 


5, R*p«h Dou 

October 1986 


P • r l©»*n »c*j> Oifanii 


7 , AulKo«( i] 


A. Kalandiya 



Pdloiminp Oifoniiclion Ko*n« *rt/ Add*m%% 

The Corporate Word, inc. 

1102 Arrott Bldg. 

Pittsburgh, PA 15222 


\ 2, Spemeiinf Apncy K#m« m*>d Add 

National Aeronautics and Space Administration 
Washington, DC 20546 


on R • por 1 No. 


10, WoiV Unil Ho, 


11, Conltoc* •• C«ont He. 

NASW-4006 

13, 7 ypm of R«perl tn^ P • • » e d Ce** •• d 

Translation 


W. Aftncy Ce^« 


li. Svpp i« m« n* or y Koi«* Translation of "K kontaktnym zadachim teorii uprugosti , " IN: 
Prikladnaya matematika i mekhanika, Vol. 21, No. 3, 1957, Moscow, pp. 389-398. 


U. Akllrocl 

Reviews certain contact problems in the two-dimensional theory 
of elasticity when round bodies touch without friction along 
most of their boundary and, therefore, Herz' hypothesis on the 
smallness of the contact area cannot be used. Derives 
fundamental equations coinciding externally with the equation in 
the theory of a finite-span wing with unknown parameter. These 
equations are solved using Multhopp's well-known technique, and 
numerical calculations are performed in specific examples. 


ORIGINAL PAGE IS 
OF POOR. QUALITY 


17. Ktjf (Selected by Author(i)} 


It* plilriUlUft 

Unlimited! 


U. S.cvfity CU»*5f. (ml tki. «.p»«0 

90. J.cvr tlf Claxlt. (•( 1U| pap.} 

2J* K«, tf P#f#« 

* Unclassified 

Unclassified 

20 
















ON CONTACT PROBLEMS OF ELASTICITY THEORY 


A. I. Kalandiya 

This article reviews certain contact problems in the two- /389 * 

dimensional theory of elasticity when round bodies touch without 
friction along most of their boundary and, therefore, Herz' 
hypothesis on the smallness of the contact area cannot be used. 

It derives fundamental equations coinciding externally with the 
equation in the theory of a finite-span wing with unknown 
parameter. These equations are solved using Multhopp's 
well-known technique^ - , and numerical calculations are 
performed in specific examples. 

1. Let a rigid die with flat symmetrical base be pr.e*»ed 
by a force acting along the die's axis into an elastic medium 
which is an infinite plane with a round opening. 

It is assumed that the die can move only forward and, in 
addition, that there are no stresses or rotation at infinity. 

The shape of the die's base (close to the contour of the 
opening) and the main vector of external forces squeezing the 
die toward the medium's boundary are given. 

The stressed state of the elastic body is sought. 

Let this elastic medium occupy the plane of variable 
C = £ + n). from which a circle with center at point £=0 with 
radius 1 is removed. We will assume that a single force of 
magnitude P, directed opposite axis , acts on the die. 


♦Numbers in the margin indicate pagination of the foreign text. 
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Boundary conditions for the task are written as (cf. [2], 


p. 429) : 


N = 0, T = 0 na L i; 


T = 0, v t = g(a) Ha La 


( 1 . 1 ) 


whereby the second of these conditions are fulfilled at contact 
arc L 2 , not given beforehand, while the first is fulfilled at 
the rest of the circumference L, L=L 1 +L 2> Here, N and T 
are, respectively, the normal and tangential components of 
external stress acting on contour L, v^ is normal (elastic) 
shift; g (cf ) , given at L 2 , is a true function of point £f=e la , 
which describes the shape of the die's base and, by virtue of 
symmetry, satisfies the condition g (tS) -g (-if) . 


Henceforth, we will assume that g(cf) has a second 
derivative along the arc of the contour which satisfied 
Gel'der's equation. 1 


To solve the problem, we will introduce Kolosov- 
Muskhelishvili ' s functions $(£) and V (£) . As we know, the 
following relation exists at the boundary of the area ([2], p. 
335). 
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4. ( 3 ) _}_ d> ( 0 ) _ 3 (t' ( 3 ) _ gsty (o) = N — iT 


( 1 . 2 ) 


According to the condition accepted above, at infinity we 
will have the following if j£| is large 


d» o = + O f- 2 ), 1(C)=f + 0 (C- 2 ) 


(1.3) 


whereby coefficients A and b are expressed as follows by 
components of the main vector (0, -P) of external forces; 


iP 


’2-(l + %) 


b = 


ixP 


2“(l+x) 


(*■ = 3 — 4v) 


(1.4) 


Here V is Poisson's coefficient. 
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The condition T=0 at L by virtue of (1.2) is presented as 


s®'(a) + o 2 «-(o) — W'(a) — (a) = 0 


(1.5) 


Hence, performing the operation 

1 ( ds 

2 J T=T 

L 

and taking (1.3) into account, we immediately find 


W(\) + (C) = b C - ;-*) + A at 

ji = Iim?[*'(C) — K _l ] at 


i:i> i 


( 1 . 6 ) 


To determine constant A, we will multiply (1.2) times 


-1 


0 d and integrate the resulting equality along L. We will 
have 

1 f A' (c) dz 




(1.7) 


On the basis of (1.6), equality (1.2) takes the form: 


<!>(;) -f <1>( 3 )_ 6(3 — 5) — A—N at L 


Hence, as above, we find 




1 r N (c) ih b 

2=rJ a — : ■: 


at i:i>« 


( 1 . 8 ) 


Differentiating (1.8) in terms of £ and taking into account 
2 

the continuity of normal stress n(tf) at boundary L, we obtain 

<1 T0=-4j’ i r§r + 4- 


(1.9) 


To determine N(o') at L, we will use the well-known formula 
((2], p. 135) 


2;ii7 f = Re Jo (ns (3) — a-i’(-) — 4(c )] \ 


( 1 . 10 ) 


5 
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If we take (1.5) into account, we obtain the following 
equation from this equality by differentiating in terms of ot 

3 = e'\ 

2* [*V+ = Rc «(* i)*!* (3) — (*+1) 3<I»'(3)-i-5«l» (3)4- 3 2 »r ( 3)J Hi, L, (1.11) 

Introducing here the limit values for the functions , 

0 *(£ > ) given by equations (1.8) and (1.9) and the expression for 
the left side of (1.6), on the basis of the latter, from 
condition (1.1) we obtain the equation 


/.— 1 


A' (z) d- 


-W[«w+^] 


xPi 

~ (1 + *) 


(°u — 3 u) = 


( 1 . 12 ) 


To determine the integration line in this singular integral- 
differential equation, we have one more relation 

(*(,)* = -/> (1 . 13) 


After finding N(c) and L 2 , we will determine stress 
functions with (1.8) and (1.6). In the particular case when the 
rigid die is a round washer inserted into an opening of the same 
radius, we will have g(<T)=0 at L 2 , and equation (1.2) will 
take the form: 


X— 1 
*+ 1 



A’(a)ih 1 { X{z)ch 



xP 

— ( 1 +*) 


K— = 0 ) 


(1.4) 


Note. Equation (1.12) (or, more precisely, the equation 
derived from (1.12) when its right side is replaced with a 
certain approximate expression) is found in V. V. Panasyuk’s 
article [3] . 

This equation is derived somewhat differently on the basis 
of the results of I. N. Kartsivadze (e.g. [2], paragraphs 125, 
126) . The author apparently was not familiar with relatively 


i 


6 



new analytical tools recently developed to solve this integral- 
differential equation (cf . below, paragraph 3) . 

Note also that the case of a round die (of the same radius 
as the opening) was considered by M. P. Sheremet'yev [4] , who 
proposed a method for constructing equation (1.14) somewhat 
different than that used here. 


However, the author's reasoning contains an omission 
because of which the equation derived in [4] is incorrect. 


As we know, in the two-dimensional deformed case, 
deformation e ^ and stresses <Tp, CT*. in polar coordinates are, 
according to Hook's law, related as follows 


l 

— /.• ( z i 



(1.15) 


where E is the modulus of elasticity. A similar formula 
suitable for the case of a generalized two-dimensional stressed 
state, takes the form: 


c -,1 — e > (:* — '<•',) 


(1.16) 


M. P. Sheremet'yev [4] uses only equation (1.16) ([4], p. 

439, equation (1.10)), in view of which his reasoning is sound 
only in the second of the two-dimensional cases. In this case, 
the author's equation [p. 441, (2.15)] indeed coincides with 
(1.14) (only this time one must take <*= (3-v)/(l+V) . One must 
take into account that, in reference [5] , force P is directed 
along axis 4* 


This note pertains also to the author's reasoning in other 
cases of the equilibrium of an infinite plane which were 
considered in the same work [4] . 
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2. Let a round elastic washer, generally with other /392 

elastic properties, be inserted into a round opening in an 
infinite medium, the opening and washer having the same 
radius. It is assumed that the washer is squeezed toward the 
surrounding material by a concentrated force applied to its 
center. (As before, we will assume that stresses and rotation 
disappear at infinity.) 

All elements relating to the washer (elastic constants, 
stress functions, etc.) are symbolized by 0, and we will write 
contour conditions of the task in the form (cf. [2], p. 207) 

N = 0, T =0 aa Iq, T — 0, tP = v 

F f (2 1 ) 

,V 0 — iT 0 = /V — iT aa L ’ 

The stress functions and y (£) , which correspond to an 

infinite plate, will, as before, be expressed by equations- 
(1.8) and (1.6) . 

Functions (£) andy o (£), which correspond to an 
elastic washer, will obviously take the form: 


% <C) = + <VQ. V, C) = 4^- + «F*(C) 

( fl o“' 215(1+*#) ’ b * = a o*» ) 

whereby functions $q*(£) and Yq*(^) ^^e holomorphs in a 
unitary circle. 

From the appropriate conditions (2.1) , we obtain the 
following, which is entirely similar to the previous 


W’.W+w.ra ici<i 

Kl<i 


where constant A is given by equation (1.7); it is assumed that 


8 


Im^Q* (0) =0. 


From (2.4), using, as before, the condition 


N (a,)=JV(_ o^) = 0 


(2.5) 


“ C* ar « end points of arc L 2 , we find 


W=4rf4^+°.( 2 -?) 


( 2 . 6 ) 


Substituting the left side for v° and v % and expressing 


the equality 


2,,0 


V ° 




f =r 


cftVc, 


P 1 da 2 


we have 


Ho Re {(•/. - 1) <I> (a) - (x -f 1) =S»' (a) + a*' (a) -f o 2 «F (=){ = 
= H Re { <x 0 - 1) cl» 4 (a) _ (x 0 + 1) ( 3 ) -f a4.; (a) + o*F 0 (a)} 


(2.7) 


Introducing the corresponding expressions given by the 
previous formulas, after certain simplification we obtain 


“w+l-l 

u 


N'(3) da 

3 — S 0 


P f A ? (3) di _ q 

2 -i J s 2 =i 


^ ( c o — ° 0 ) 


( 2 . 8 ) 


where 

/C— : 1 (1 ~ 2V) (1 + V) E ° ~ (1 ~ 2V(>) (1 + V ° } £ 

2 (1 - v2) /T 0 + (1 — vj) £ 

(1-v 2 )£q j - 1 (1+v)^0 + (l+VQ)g 

P (i-v2)£ 0 + (i-v 2 )£ ’ * 2 (! _ v 2) £ 0 -f- (1 — vjj) E 
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(2.9) 

(2.10) 


whereby E represents the modulus of elasticity. The contact arc 
I »2 is, as before, unknown, by virtue of which relation (1.13) 
should be attached to equation (2.8) as before. 


When the washer and the plate are made of the same material 
(Vq=Vp Eq=E) , coefficient k in (2.8) becomes zero, as a 
result of which this equation is solved in closed form. The 
solution (closed) to the problem in this case was found in [5] . 
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Moving in equation (2.6) to the limit at E^oO, we come to 
the case of an absolutely rigid washer -- equation (1.14). 


If E=<», we have the case of a round, elastic washer 
inserted into an opening of equal radius in an absolutely rigid 
plate pressed toward it by a concentrated force applied to the 
center of the washer. Equation (2.8) gives 


x— l 
x+l 




A" (c) rfo 


P 

1+x) 


(°o — 3 0 ) 


whereby X. is the washer's elastic constant. 


( 2 . 11 ) 


Equation (2.8) is an example of an equation for a two- 

dimensional contact problem covered without any limiting 
3 

assumptions. 

3. Let us add to equation (2.8) a new variable 
introduced by the relation 


. x — i'i a COS a, 

3 1 x + if ’ " i+sina*’ 


o* = e’*» 


(3.1) 


The previous function, as we know, gives the conversion of 

the unitary circumference l<r| =1 into real axis x, transferring 

arc L 0 to segment [-1, 1] ; the point -O'*, O’* on the 
^ 4 

circumference is transferred respectively to points x=-l, x=l. 

After the obvious transformations, we will have 

K/ t \ 1 f 1 A” (t)dt P p f N (t) dt » xt-p 

2 *J t — x = (**+?*) J P + p ^(*2+12)2 

N(l) = N( — 1)=0 

whereby, for simplicity, we again assume N(tf)=N(x). The 
equality (1.3) is transformed to 

i 

r J2 fit 

2,3 J N{t)it = P 

—1 


(3.2) 


(3.3) 


1 


10 


/394 


If we disregard unknown parameter p, equation (3.2) is 
the familiar equation of the theory of a finite-span wing, which 
is usually written as: 


r(g) 

B(x) 



r (t) dt 

t — X 


=/(*) 


(3.4) 


Here B and f are assigned functions at segment [-1, 1] , 
wherein B(x) is zero nowhere except the ends of this segment, 
while r (x) is the desired function subordinate to condition 

r (l) = f(-l)=0. 


Many articles have been devoted to the wing theory equation 
(cf. e.g. [7, 8]). Among the countless methods of numerical 
solution devoted to this equation, the most successful from the 
standpoint of practical applications is Multhopp's (direct) 
method [1], which, because of its simplicity and efficiency, is 
still considered the best mathematical device for aerodynamic 
calculation of a wing. 

According to those methods, the desired T takes the form of 
a trigonometric interpolation polynomial 


L[I\ x]=S(-ir 1 r(x,) 


s+i n \ sin (n -f- 1) a sin »» 


Ie=1 


COS i) — cos 


y. r (x k ) y sin sin nt!> 
n + 1 , 


k=l 


m=zl 


(3.5) 


; = cos&, x* = cos&n., 8f = n ^ (0 < 9 ^ c; fc — 1.. 


n) 


and, after application of a certain squaring formula to a 
special integral on the left side of (3.4), this equation is 
replaced by a system of linear equations of relatively 
approximate values of the desired function in the given 
(Chebyshev) nodes x^. This system takes the form 


+ &».) I\. =/„ + ^ b*T k (m = 1 «) 


(3.6) 
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where B +B(x ), f =f(x ); b . is also known, so that 
m m m m ink 

b m ^=0 at m-k =2, 4... Substituting the solution to this 
system into the right side of (3.5) instead of f"^, we obtain 
the approximate solution to the equation (3.4). 


Multhopp's article [1] proved that consecutive 
approximation method for (3.6) always converges if B(x) is 
negative [1]. With negative B(x), which occurs (as can be seen 
from equation (3.2)) in applications to elasticity theory, 
iteration method as applied to system (3.6) also yields a 

g 

converging process, provided that B(x) satisfied the condition 


max 

*€(— 1,1 ] 




< ~ 


(3.7) 


We will apply this method to solving the equations derived /395 

here. 


Solving (3.2) and (3.3) together, we will, according to 
(3.5) , find the approximate solution to (3.2) in the form 


n (x) = 2 V ~r ~ 2 (*). a - = sin w,) - 


TT , . sin mft 
<7 --> 


(3.8) 


The function U , (x) is a polynomial of the x-th order of 
m-i 

m-1 (a second-order Chebyshev polynomial) . 

To calculate the desired values for the problem, it is 
convenient to shift from the physical plane of variable ^ to the 
plane of variable z=x+iy with the relation 

<=‘St (?=t rfc) < 3 - 9 > 

which accomplishes the conformal transformation of the circle 
\^|<1 to the upper half-plane IM z^O. After transformations, we 
obtain the formula 


* 
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(3.10) 


i i 

\ [N(z)d<s 1 f N (x) dx . 3 CN (x) dx 

2~i J 5 — r — 2-i J a: — z + 2 r J *2 + ?* 

-l -l 

But the Koche integral on the right side of (3.10) is 
calculated in final form if N(x) takes the form of (3.8) . 
Consequently/ after Multhopp's method is used to solve the 
equation in the contact problem, the corresponding functions for 
stress are determined in closed form. 


Example 1. An absolutely rigid round washer squeezed 
against an elastic plate (external problem). The problem's 
equation is derived from (1.1) by substituting (3.1) or the 
limit transition in (3.2) at Eq*«o. We will write this 
equation as 

(3.11) 


X— 1 


-13 if 

Hr 2 — 3‘- iV ( r ) 2s .1 
-1 


K’ ( t ) dt 


T-q 


2x3 P 


•A _ 3 * 


x- -f 1 « ( 1 -i- x ) (** + J 1 *) 2 


where 

1 f N ( t ) dt 

9— - J 12-1-32 

—l 


(3.12) 


Let us remember that parameter introduced into (3.11) 
and characterizing the size of the contact area, is unknown and, 
in addition, the right side of the equation contains the still 
undefined constant q in the form of functional (3.12). In the 
rest, (3.11) coincides with (3.4), so that BJx) } 0 at (-1, 1], 


Equality (3.3) is rewritten as: 


23 f / 2-32 
c = P J (<2 4- 32)2 A 


((> _2_ J}2)2 A (0 dl — 1 


(3.13) 


In our (symmetric) case of equation (3.11), if n is odd, 
system (3.6) will take the form: 

— l 




+ 1 cos 2 o M 4- y, 


3 2 + 


»~| A in 


v ' ?2<l 1_ cos2 ~ 

— z J — co S 2 S (1 + x) (cos 2 H*-}-.* 2 ) 


(3.14) 


»=1 


13 


where 


I'ir — r„ +1 _i- 
B,Hh — bulk -)- b„ i, H4-1 — t 


(m— 1,2 1/2 ("+ 1 )) 

lt= 1, . . . , '/;<» — 1). bm,'UU>+h 


(3.15) 



Solving this system for a certain parameter {3 (e.g. 
, we will find unknwon in the form 


at 
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(* = 


(3.16) 


whereby will be known. The corresponding 

approximate solution to N(x) will obviously contain constant q, 
defined as a result of (3.12). Then, after N(x) is calculated, 
the left side of equality (3.13) is determined. To determine q 
and c(^), a Gaussian squaring formula is used (e.g. (9], p. 617). 

1_ vi -sin 8 t A't 

9 « -f 1 Zj cos- fti. -J- 3* 

t=i * ■ ' 

c (3) — 2 ^~ V sin 8 * (cos 2 ~ . 32 ) Nk (3.17) 

n + 1 ^ (cos* + £*)* P ~ 1 

The value of ^ which we have taken and the N(x) 
correspnding to it will generally not satisfy equality (3.13). 
Therefore, selecting new $ values, we will repeat our 
calculations until this equality is satisfied with necessary 
accuracy. As a result, we will have certain approximate values 
for ^ and N(x) for a given n. Then other n's are taken and 
these calculations repeated. 

In system (3.14), we will first assume n=7 and will set 
Poisson's coefficient equal to 1/3 for further calculations. In 
expanded form the resulting system appears as:^ 

£4 (cos* #, + £2) + 5 - 2262 ] = 

_ _ 5 a cos*»i— £* P 

= 1.9142 (V, + 0.1464A 4 + cos2&i _j_ + 4 ^ cos 2a,+p2)2 „ 
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—0.9142 W 2 -j- 0.8536 + cos2 ^ _|1 p2 + 4 ? ( C os2»j+?2)2 it 

[ 4(cos4 2 -t-^) + 2 - 8284 ]^ = 

_ 032 , 5 cos2»~ — 32 p 

- 1 0360 ^1 + 1- 1944 N 3 + cos2 y 2 _j_ 3-2 + 4 3 ( COS 2B 2 +32)2 T 

(_ 4 (cos 2 ft 4 4 - 3 2 ) + 2 000 °] N i = 

- - . «32 5 COS2&4 — 32 p 

- 0.1121 JVj-f 1.5774 yv 3 + cos2 ^ + ' j2 + 4 3 (cos2 » 4 4 - 32)2 - 

®jt == "g“ 1 A - = 1 , 2, 3, 4 


(3.18) 


Solving this system together with (3.17), we will have 
(values for q and N k [1] in (3.16) as well as N k are given 
in P/t\ ) 8 


3 = 1.20886, q= -0.52422 ( 3.19) 

= 0.48035 q — 0.29067 = —0.54248, A* = 0.91050 q — 0.63176 = —1.10906 
N z — 1.23317 q— 1.00926 = — 1.65571, A'j = 1.35900 j— 1.20474 = —1.91715 ( 3.20) 

The value of the following polar angle corresponds to the /397 
^ which has been found 

= -10"48’12" (3.21) 

with which contact arc is calculated. Given these values 
for the unknowns, the left side of (3.17) equals 

c (3) =0.99947 

9 

For maximum pressure in our approximation, according to 
(3.20) we will have 


' V Lo^'V 4 = —0.61025 P 


(3.22) 


Later in (3.14) we will use n=15. This system of eight 



equations (with five separate unknowns in each) will be solved 
by consecutive approximation method, wherein we will use the 
values in (3.20) with unknown q as the zero approximation of 

N 2' N 4' N 6' N 8* 

Solving this system with the $ found together with (3.17), 
we obtain Multhopp's second approximation 

«7 = — 0.52418 ( 3.23) 

.Vj = 0.24397 q — 0.14220= — 0.27008, N 2 = 0.48054 g — 0.29080 = — 0.54269 

jV 3 = 0.70465 q— 0.45269 = —0.82206, A T 4 = 0.91061 g— 0.63180 = —1.10913 

yV 5 = 1 .09038 q — 0.82382 = — 1.39537, A c = 1.23312 g— 1.00927 = -1.65565 

jV 7 = 1 .32624 g - 1.15058 = -1.84577, N 8 = 1.35877 g— 1.20437 = — 1.91663 (3.24) 


For these" values of Q and N k (k=l, 
of (3.17) gives 


C(’1) = 0.99935 


8) , the left side 


For maximum stress, we have 


N 1^=0 ~A^g = — 0.61008 .P (3.25) 

As we can see, the desired values in the first and second 
approximations differ little from one another. 

2. A round, elastic washer, squeezed against an opening 
in a rigid plate (internal problem). The appropriate equation 
is derived from (2.11), if we introduce (3.1) into it. It takes 
the form: 

* — 1 ? „ , 1 V N’ (t) dt 2.3 P X*- — P 

x + lxS + .32"W-2eJ t-x -=(l+x)(x2+^)2 (3.26) 

This equation (with the stipulation made above relative to 
^) also coincides with (3.4), so that, this time, B(x) 0 at [-1, 
1 ]. 
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System (3.6) for equation (3.26) externally differs little 
from system (3.14). It takes the form: 


t ~ 1 

7. -f- 1 cos- »■ -r -r ' 


bmm 






•Snii'Yjt -f r 


25/> 


cos 2 8* — a 2 


F=] 


( 1 + x) (cos 2 8* + > 2 ) 2 


("»=1 Vs (»+!)) 


(3.27) 


This system should be solved, as before, by satisfying 
equation (3.17) and simultaneously determining parameter At 
n=7, we will have the following values for the desired parameters 
(N k is given in P/H) : 


3 = 0.91612 (3.28) 

JVj = — 0.37480, iV 2 = —0.85751 , A 3 = — 1.47396, 7V 4 = -1.84453 ( 3.29) 

Given these values, c(^)=l. 00034. 

From this we obtain the following for the desired value of /398 
angle or, which corresponds to the end of the contact arc 


a, = 5° 00' 46" 


(3.30) 


and maximum pressure will equal 

/V I *=o~iV t = — 0.58713P (3.31) 

In (3.27) we will then assume n=15. Solving this system 
with (3.28), we obtain^ 

= —0.17990, .V 2 =— 0.37458, iV s = —1.15746, ;V 8 = — 1.47378 
N 3 = - 0.59686. .V 4 = — 0.85727, ;V 7 = — 1.73637. N s == —1.84215 (3.32) 

Theft c(^) =0.99905. Specifically, for maximum normal stress 
in the second approximation we will have 

N ! *=o~.V= —0.58638 A 
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we can, as 


According to the values found for N(x) and 
indicated above, calculate all other desired values for the 
corresponding contact problem. Specifically, it is easy to 
calculate normal shifts which will satisfy the contour condition 
of the problem with a certain error.^ 

In conclusion, let me note that these calculations were 
done at the USSR Academy of Sciences' Computing Center. I will 
take the opportunity to thank the directorate of the Computing 
Center for its assistance, as well as Ye. S. Bogomolova and T. 

M. Kopylova of the Computing Department for their rapid and 
accurate work. 
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FOOTNOTES 


1. The last condition of (1.1) should have been written more 
accurately: Vp=g(cf)+c sin cx., where c is the forward 
movement of the die; but one can do without this shift 
since it can be eliminated by the rigid forward shift of 
the entire system. 

2. It is assumed that the die has no angular points at the 
contact boundary. On this assumption N(cj) will be a 
continuous function at L, becoming zero at the ends of the 
contact arc. 

3. Contact problems of this type, which result in equations of 

the same structure, were first considered apparently in I. 
Ya. Shtayerman's work (e.g. [8]). In this article, we 

direct our attention primarily to the possibility of 
efficiently solving these problems. 

4. <T* denotes the end point of arc L 2 for which Recf*>0. 

5. The functions B(x) are exclusively such in application to 
wing theory. 

6. Justification of Multhopp's method and its application in 
problems discussed here are found in the author's 
dissertation ("Certain Mixed Problems of Elasticity 
Theory," Steklov Mathematical Institute, USSR Academy of 
Sciences, 1955) . 

7. Reference [1] has tables for calculated values of b mk at 
n=7.15.31. 

8. Here and henceforth we will require that relative error in 
(3.17) not exceed 0.1%. 

9. The solution to N(x) which corresponds to n=7 is assumed, 
according to Multhopp, to be the first approximation. 

10. Given our and the Poisson's coefficient taken above, 
function B(x) in (3.26) satisfies condition (3.7). 

However, we must state that consecutive approximation 
method converges much more slowly in this case than in the 
previous case. For this reason, the method cited, when 
applied to system (3.6) given negative B(x) {even with 
(3.7)], may not always be preferable. 

11. ‘For example, the absolute error for/i/P(Vp+ 

+d 2 v p /da 2 ) for the external problem does not exceed 
0.0072. 


I 
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